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Abstract 

The closely related protozoan parasites Toxoplasma gondii and Neospora caninum display similar life cycles, subcellular 
ultrastructure, invasion mechanisms, metabolic pathways, and genome organization, but differ in their host range and 
disease pathogenesis. Type II (y) interferon has long been known to be the major mediator of innate and adaptive immunity 
to Toxoplasma infection, but genome-wide expression profiling of infected host cells indicates that Neospora is a potent 
activator of the type I (a/[5) interferon pathways typically associated with antiviral responses. Infection of macrophages from 
mice with targeted deletions in various innate sensing genes demonstrates that host responses to Neospora are dependent 
on the toll-like receptor Tlr3 and the adapter protein Trif. Consistent with this observation, RNA from Neospora elicits TLR3- 
dependent type I interferon responses when targeted to the host endo-lysosomal system. Although live Toxoplasma fail to 
induce type I interferon, heat-killed parasites do trigger this response, albeit much weaker than Neospora, and co-infection 
studies reveal that T. gondii actively suppresses the production of type I interferon. These findings reveal that eukaryotic 
pathogens can be potent inducers of type I interferon and that related parasite species interact with this pathway in distinct 
ways. 
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Introduction 

The eukaryotic phylum apicomplexa is comprised of over 5000 
species of parasitic protozoa that infect a wide range of animal 
hosts and cause significant disease in both healthy and immune- 
compromised individuals. This phylum includes Plasmodium, the 
causative agent of malaria; Cryptosporidium, recently recognized to 
be a leading cause of pediatric diarrheal disease in the developing 
world [1]; and Toxoplasma, a ubiquitous parasite that causes 
potentially fatal congenital disease and represents the second 
leading cause of food borne mortality in the USA [2]. Other 
apicomplexa, including Babesia, Neospora and Theileria are impor- 
tant pathogens of cattle, while Eimeria species are a leading concern 
in the poultry industry. Despite being important causes of human 
and animal disease worldwide, there is relatively little information 
about how these pathogens are recognized by the host innate 
immune system following invasion of target cells, nor is it known 
whether infections with distinct protozoa trigger identical, partially 
overlapping or completely distinct innate immune responses. 

The profound susceptibility of IFN-y-deficient mice to a wide 
range of protozoan parasites, including Toxoplasma gondii [3,4], 
Cryptosporidium parvum [5] , Leishmania major [6] and Trypanosoma cruzi 
[7], has given rise to the widely-held view that type II interferon 



(IFN-y) is responsible for controlling protozoan infections, in 
contrast to type I interferons (IFN-ot/P) which are primarily 
associated with control of viral infections. As a consequence, the 
vast majority of studies on innate immunity to protozoan parasites 
have focused on IFN-y produced during acute infection by NK 
cells [8,9], and interleukin- 1 2 produced by macrophages [10], 
neutrophils [11,12] and dendritic cells [13]. It is unclear whether 
this accurately captures the complete picture of innate immune 
pathways triggered by infection with protozoa. Moreover, while 
tremendous progress has been made in identifying viral and 
bacterial ligands recognized by innate pattern recognition 
receptors [14], relatively few studies have explored the mecha- 
nisms and consequences of innate recognition of intracellular 
eukaryotic microbes [15], in part because the conserved signatures 
of bacteria and viruses are not generally thought to be present in 
eukaryotic cells. 

We have employed a comparative genomic approach to 
investigate innate immune signaling triggered by two closely 
related but distinct apicomplexan species, Toxoplasma gondii and 
.Neospora caninum, both of which are obligate intracellular parasites 
that initiate infection in the gastrointestinal tract of animals. 
Neospora was only recently recognized as a distinct protozoan 
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species [16], and although it is thought to have diverged from 
Toxoplasma ~28 million years ago, the two species share nearly 
indistinguishable ultrastructure [17], metabolism, and 1:1 syntenic 
orthologs for >90% of the parasite genome [18]. Despite these 
similarities, important biological differences distinguish the two 
species with respect to virulence factors, host range, and 
pathogenesis. Toxoplasma undergoes sexual reproduction in felidae 
(cats) and asexual reproduction in virtually any mammalian host, 
while the Neospora sexual cycle takes place in canidae (dogs), where 
it demonstrates a propensity for vertical transmission from mother 
to fetus [19], and infects a more limited range of intermediate 
hosts [20,21]. Genomic and transcriptomic analysis of these two 
parasites reveals that Neospora and Toxoplasma have evolved distinct 
repertoires of surface antigens and secreted kinases since 
divergence from a common ancestor [18], suggesting the ability 
modulate or interact with host cells in unique ways. We capitalize 
on the close genetic and ecological relationship between these two 
parasites to carry out comparative profiling of the host responses to 
infection. Our results demonstrate that Neospora is a potent inducer 
of innate IFN-oc/ (3 responses, and that Toxoplasma has evolved the 
capacity to suppress this response. 

Results 

Differential induction of innate immune signaling by 
genetically similar parasite species 

In order to compare host transcriptional responses during 
Neospora and Toxoplasma infection, whole-genome expression 
profiling was carried out using human fibroblasts infected with 
either Neospora caninum (NcLiv strain) or Toxoplasma gondii (GT1, 
Prugniaud, or VEG strains, selected as representatives of the three 
dominant genotypes observed in North America [22,23]). 822 host 
genes were identified as differentially transcribed (^2-fold, 
FDR<5%) during infection with any of these four parasites 
relative to uninfected host cells (Fig. la and Table SI). Grouping 
parasites based on the host response indicated that Toxoplasma 
strains known to exhibit low virulence in mice (PRU, VEG) elicit 
host response profiles more similar to each other than to the 
virulent GT1 strain. Neospora induced a host response distinct from 
all strains of Toxoplasma, but more closely resembling low-virulence 
Toxoplasma strains, indicating that these two parasite species 
perturb the host transcriptome in distinct ways. 

Hierarchical clustering of the 822 genes revealed a distinctive 
set of 66 genes induced only during Neospora infection (Fig. la, 
asterisk; Fig. lb), and Gene Ontology (GO) analysis [24] indicated 
that these are significantly associated with antiviral responses 
(Fig. lc). Antiviral genes specifically induced by Neospora infection 
included interferon-JM (Ifnbl) as well as several key regulators of 
type I interferon responses (Fig. lb, arrows) [25-27], including the 
transcription factor IrfJ (a master regulator of the antiviral 
program), Oasl and Oas2 (members of the 2'-5'-oligoadenylate 
synthetase family that activate RNase L to initiate degradation of 
cellular RNA as an innate defense mechanism during viral 
infection), and Mxl (a GTP-binding protein that is essential for 
antiviral responses to influenza infection). To test whether the 
transcriptional program induced by Neospora is functionally 
antiviral, uninfected human fibroblasts or cells infected with 
Neospora or Toxoplasma were challenged with GFP-tagged Vesicular 
Stomatitis Virus (VSV), a standard virus used in bioassays for type 
I interferon [28] . Imaging these cells demonstrated that while VSV 
grew normally in naive and Toxoplasma-infected cells, the virus was 
dramatically restricted in jV"««/>ora-infected cultures (Fig. Id). 

Viruses are known to infect various parasitic protozoa [29], in 
some cases influencing host transcriptional responses and promot- 



ing disease [30]. Several arguments indicate thatindicate t a 
cryptic Neospora virus is unlikely to be the source of the type I 
interferon response noted above. Analysis of raw trace reads from 
the Toxoplasma and Neospora genome projects [18] failed to identify 
evidence of viral DNA. Similarly, deep sequencing of both 
polyadenylated mRNA and small non-coding RNA libraries 
prepared from Neospora and Toxoplasma infected human cells 
yielded no significant evidence of viral RNA, based on mapping to 
NCBI RefSeq virus genomes and comparison with uninfected host 
cells (data not shown). Finally, total RNA isolated from Toxoplasma 
or Neospora and treated with S 1 nuclease to degrade single stranded 
nucleic acid showed no evidence of double stranded RNA species 
smaller than 1 2 kb after staining with ethidium bromide (Fig. S 1 a). 
In both parasite species, SI nuclease-resistant bands greater than 
12kb (Fig. SI a) were DNase sensitive (Fig. Sib). In contrast, this 
approach readily detected double stranded RNA from the 
endogenous RNA virus, LRV1, present in Leishmania guyanensis 
(Fig. SI, arrows) [31]. Taken together, these data suggest that a 
parasite molecule present in Neospora, rather than an endogenous 
virus, is responsible for induction of a type I interferon response. 

Toxoplasma actively suppresses host cell induction of a 
type I interferon response 

Using Mxl transcript abundance as a read-out for host response 
to parasite infection, three strains of Neospora and four strains of 
Toxoplasma were tested for induction of the type I interferon 
program. All Neospora isolates were potent inducers of Mxl 
expression, while no strain of Toxoplasma tested induced this 
response (Fig. 2a). Induction of Mxl correlated with parasite 
inoculum (Fig. S2a), suggesting that infected cells (rather than 
uninfected bystanders) are the likely the source of this response. 
Because N. caninum is a common pathogen of cattle rather than 
humans [20,21], primary bovine fibroblasts were also infected with 
Neospora, demonstrating that both human and bovine cells 
upregulated Mxl expression following infection (Fig. S2b). 

Many responses to Toxoplasma infection require active invasion 
of the host cell, including activation of the transcription factors 
STAT3 and EGR2, and suppression of proinflammatory signaling 
[32-34], indicating that these phenotypes are mediated by parasite 
factors secreted at the time of invasion or during intracellular 
replication. To test whether induction of the type I interferon 
profile also requires active infection, parasites were heat-killed 
prior to inoculation. Heat-killed Neospora retained their ability 
induce an antiviral program, even exhibiting enhanced Mxl 
induction compared to live parasites (Fig. 2b). Surprisingly, heat- 
killed Toxoplasma strains also induced Mxl expression, albeit at a 
level that averaged 1 1 -fold weaker induction than seen with heat- 
killed Neospora strains. These data demonstrate that although both 
parasite species harbor a factor able to trigger the type I interferon 
response, Neospora elicits a much more robust response. 

Given that innate induction of type I interferon is usually 
triggered by Toll-like receptors or RIG-I-like receptors (RIG-I and 
MDA5) that are located in the endo-lysosomal pathway or free in 
the cytosol, respectively, we were surprised that killed parasites not 
only induced a type I interferon profile, but did so more strongly 
than live parasites. We hypothesized that heat-killed parasites must 
get internalized by fibroblasts in order to access innate receptors 
that elicit this response. To test this idea, we used confocal 
fluorescence microscopy to visualize heat-killed Toxoplasma ex- 
pressing mCherry in cells labeled with LysoTracker to identify 
lysosomes (Fig. 2c). While live Toxoplasma could be seen replicating 
in LysoTracker-negative vacuoles 18 hours after infection, heat- 
killed parasites were internalized and localized to LysoTracker 
positive compartments in the cell (Fig. 2c, arrows). 
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Figure 1. Differential induction of innate immune signaling by genetically related parasite species, a) Microarray-based expression 
profiling of human foreskin fibroblasts (HFF) infected with Neospora caninum (NcLiv isolate) or Toxoplasma (GT1, Prugniaud or VEG strains). Heat map 
shows hierarchical clustering analysis of 822 genes differentially regulated relative to uninfected cells by >2-fold (FDR<5%), in any of these 
experiments. Each row in heatmap represents average of duplicate (NcLiv and GT1) or triplicate arrays (VEG, PRU, uninf). Color pattern on heatmaps 
represents column Z-score. b) A cluster of 66 genes (Fig. 1a, asterisk) induced by Neospora but not any strain of Toxoplasma, including several well- 
known type I interferon response genes (arrows), c) Gene Ontology (GO) enrichment analysis of 66 Neospora-induced genes. Bar graph shows fold 
enrichment for top five GO Biological Process terms enriched at PsO.05 and represented by >5 genes. Number of genes and GO term name is shown 
at the right of each bar. d) Fluorescence micrographs of uninfected HFF cells, or cells infected with Neospora or Toxoplasma, then challenged with 
GFP-tagged vesicular stomatitis virus. Representative images are shown. Experiment was repeated three times with similar results. 
doi:1 0.1 371 /journal.pone.0088398.g001 



The observation that heat-killed, but not live, Toxoplasma induce 
Mxl expression suggested that live T. gondii actively suppressed this 
response. To test this hypothesis, human fibroblasts were first 
infected with live Toxoplasma (which does not induce Mxl 
expression), and then challenged 2 hr later with heat-killed 
Neospora (which do induce Mxl). Supernatants were collected at 
24 hr after challenge, filtered to remove parasites, and transferred 
to naive fibroblast cultures to assay for type I interferon (Fig. 3a). 
While media from Neospora treated cultures was a potent inducer of 
Mxl expression, media from cells infected with Toxoplasma prior to 
Neospora exposure failed to elicit Mxl expression (Fig. 3a). These 
results demonstrate that Toxoplasma infection suppresses innate 
immune induction of type I interferon. 

To test whether the suppressive effect observed in co-infection 
experiments was a cell-autonomous phenotype evident only in cells 
infected with Toxoplasma, conditioned media was recovered from 
Toxoplasma infected human fibroblasts, filtered to remove parasites, 
and used to pretreat naive fibroblasts (Fig. 3b). Pretreatment of 
cells did not prevent subsequent infection with Neospora (not 
shown), but blocked Neospora induction of the antiviral genes IRF7 
and HERC5, 3-fold and 4-fold, respectively (Fig. 3b). Supernatants 
were filtered again and transferred to new cells to test for ability to 
restrict a viral challenge. As expected, cultures treated with 
supernatants from Neospora infected cells were completely refrac- 
tory to a challenge with VSV-GFP (Fig. 3c, upper right). In 
contrast, supernatants recovered from cultures that were pretreat- 
ed with Toxoplasma conditioned media prior to Neospora infection 
showed impaired virus restriction (Fig. 3c, lower right). Taken 
together, these data show that Toxoplasma infected cells produce a 



soluble factor that suppresses the innate induction of type I 
interferon. 

The macrophage response to Neospora infection is 
dependent on signaling through the type I interferon 
receptor 

In order to determine the extent to which host transcriptional 
responses to Neospora and Toxoplasma (Fig. 1) are driven by type I 
interferon, bone marrow-derived murine macrophages were 
generated from either wild-type mice or mice lacking the type I 
interferon receptor subunit Ifnarl [35] . Transcriptional profiling of 
wild-type (WT) and Ifiiarl -deficient macrophages infected with 
Neospora or Toxoplasma identified 833 genes differentially regulated 
(&2-fold, FDR£5%) compared to uninfected cells (Fig. 4a, Table 
S2). Although Toxoplasma elicited a similar host response on both 
cell types, the host response to Neospora was severely abrogated by 
the loss of type I interferon signaling (Fig. 4b), demonstrating that 
signaling through this receptor constitutes the major pathway 
responsible for shaping the transcriptional response to N. caninum. 

Hierarchical clustering identified two distinct clusters of co- 
regulated genes (Fig. 4a). Cluster 1 (212 genes) was induced by 
Toxoplasma, but not Neospora, and was ijfoar/-independent. Gene 
ontology analysis showed that this cluster was strongly enriched for 
terms relating to wound healing (Fig. 4c), and arginase-1, Chi313, 
Earl 1 and I14il were amongst the most strongly induced genes in 
this cluster (Fig. 4d) - all well-known markers of alternative 
macrophage activation, which have recently been shown to be 
induced during infection with certain strains of Toxoplasma [36] . In 
contrast, cluster 2 (199 genes) was strongly induced by Neospora, yet 
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Figure 2. Active invasion is not required for innate recognition of parasites. QPCR analysis of the expression of the antiviral gene Mxl 
following (a) infection of HFF cells or (b) treatment of cells with heat-killed strains of Toxoplasma (open bars) or Neospora (shaded bars), (c) confocal 
fluorescence microscopy of HFF cells treated with live (top row) or heat killed (bottom row) 7bxop/asma-mCherry. Lysosomes are stained with 
LysoTracker dye. Representative images are shown. Error bars indicate standard deviations for three biological replicates; * = P<0.01. Experiments 
were repeated three times with similar results. 
doi:1 0.1 371 /journal. pone.0088398.g002 



only weakly induced by Toxoplasma, was Ifnarl -dependent (Fig. 4a), 
and was enriched for GO terms related to immune cell function 
and inflammation (Fig. 4e). The top ten most strongly induced 
genes in this cluster included key STAT 1 -dependent genes such as 
Cxcl9 and Cxcll 0 as well as Ir/7, a master regulator of antiviral gene 
transcription (Fig. 4f). Taken together with our data on human 
fibroblasts and bovine fibroblasts, these results show that Neospora 
induces a strong type I interferon signature in cells from diverse 
host species as well as immune and non-immune cell types, 
suggesting a conserved innate response to these parasites. 

The type I interferon response to Neospora is mediated 
by signaling through toll-like receptor 3 (TLR3) and the 
adaptor protein TRIF 

To identify the host pathway responsible for parasite-mediated 
induction of type I interferon, and downstream responses to this 



cytokine, bone marrow-derived macrophages were prepared from 
mice genetically deficient in various toll-like receptors or their 
adaptor proteins (Fig. 5). Neospora infection of Myd88 deficient cells 
still resulted in induction of Mxl expression, ruling out a role for 
TLR1, 2, 5, 6, 7, 8, and 9 [37]. TLR4 is capable of inducing type I 
interferon production via either MYD88 or the alternative adapter 
TRIF [38], but infection of Tlr2/Tlr4 double knock-out cells 
induced antiviral gene expression, ruling out a role for TLR4 as 
well. In contrast, cells lacking Tlr3, or lacking both Mjd88 and Trif, 
showed no induction of Mxl by live Neospora. Tlr3 ' and 
Myd88~ ' f ~ Trif ' ~ cells also failed to respond to heat-killed 
parasites (not shown), demonstrating that both live and dead 
parasites are sensed through this pathway. Since TRIF is required 
for TLR3 signaling [39,40], these data identify TLR3 - known to 
recognize double stranded RNA [41] - as the host receptor 
responsible for macrophage sensing of Neospora and type I 
interferon responses. 
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Figure 3. Toxoplasma actively suppresses anti-viral host responses via a soluble factor, (a) co-infection assay in HFF cells showing Mxl 
gene expression in Toxoplasma infected cultures (Tg), cultures inoculated with heat-killed Neospora (HK Nc; striped bar), or infected with Toxoplasma 
for one hr before challenge with heat-killed Neospora (black bar), (b) QPCR analysis of IRF7 (open bars) and HERC5 (striped bars) expression in HFF 
cells pre-treated with fresh media or conditioned media from T. gondii infected cells (1 ° treatment) prior to infection with N. caninum (2" treatment), 
c) HFF cells were pretreated with supernatants from experiment in panel a, then infected with VSV-GFP. Representative fluorescence micrographs are 
shown. Error bars indicate standard deviations for three biological replicates; * = P<0.01. Experiment was repeated two times with similar results. 
doi:1 0.1 371 /journal.pone.0088398.g003 



Although TLR3 is primarily localized within endosomes, some 
cell types also express this receptor on the plasma membrane 
[42,43]. To determine whether surface or endosomal TLR3 is 
responsible for sensing live and heat-killed Neospora, macrophages 
were pretreated with bafilomycin Al (Fig. 6), a specific inhibitor of 
vacuolar ATPase that blocks endosomal acidification and subse- 
quent fusion [44]. Although both live and heat-killed parasites 
induced the expression of Mxl and hf7 in macrophages treated 
with DMSO alone (Fig. 6a), this response was abolished in 
bafilomycin-treated cells, suggesting that endosomal TLR3 is 
responsible for recognizing both live and dead parasites in 
macrophages. 

Since RNA is the prototypical ligand for TLR3 [41], parasite 
RNA was tested for its ability to activate this pathway. RNA 
extracted from Neospora or uninfected human cells (to control for 
any contamination of parasite material with host cell RNA) was 
added directly to Myd8S~'~ macrophage cultures (to exclude 
activation of TLR 7, 8 and 9, which also bind nucleic acids). None 
of the RNA preparations induced Mxl or Irf7 expression, although 
the synthetic TLR3 ligand poly(LC) resulted in 15-fold induction 
(Fig. 6b; white bars). Because it is unclear how efficiently 
extracellular RNA is taken up into the endosomes where TLR3 
resides, parasite RNA was also complexed with the liposomal 
transfection reagent l,2-dioleoyl-3-trimethylammonium-propane 
(DOTAP). When complexed with DOTAP, RNA from Neospora, 
but not Toxoplasma, strongly induced both Irf7 (Fig. 6b and 6c) and 
Mxl (not shown). Thus, it appears that naked Neospora RNA is 
capable of triggering 7l'i)«/i9<?-independent type I interferon 
responses in macrophages, while RNA from Toxoplasma is not. 
Although the use of Myd88~'~ cells in these experiments helped to 
narrow the range of TLRs that could be acting as receptors for 



Neospora RNA, it does not rule out a role for cytosolic RNA sensors. 
Therefore, we tested specifically whether TLR3 was the receptor 
for transfected parasite RNA by pre-incubating Myd88 1 
macrophages with a small molecule competitive inhibitor of 
dsRNA binding to TLR3 [45]. In the presence of inhibitor, IRF7 
induction by Neospora RNA was completely abolished, demon- 
strating that this is a TLR3-dependent process (Fig. 6c). 

Discussion 

Despite the importance of protozoan parasites as causes of 
human and animal disease, relatively little information is available 
on innate signaling pathways triggered by these organisms, in 
contrast to the extensive literature on innate responses to viral and 
bacterial pathogens [14]. Glycosylphosphatidylinositol (GPI)- 
anchored surface antigens (found in many parasite species) can 
bind to TLR2 and/or TLR4 [46]; Plasmodium hemozoin can 
activate the inflammasome [47] and (when bound to DNA) TLR9 
[48,49]; and Toxoplasma profflin binds to TLR1 1 [13] and TLR 12 
[50,51]. These interactions, in turn, trigger the production of pro- 
inflammatory cytokines, including IL-1, IL-6, TNF-a and IL-12. 
To this list, we can now add the induction of type I interferon 
responses (Fig. 1 and Fig. 4) via TLR3 and TRIF (Fig. 5), and our 
in vitro data suggest that parasite RNA is a potential trigger for this 
response (Fig. 6). 

Infection by protozoan parasites has typically been associated 
with induction of type II (y) interferon, while type I (oc/P) 
interferon is prominent in viral infection, but type I responses to 
protozoa are not completely unprecedented. Toxoplasma infection 
elicits type I interferon from plasmacytoid dendritic cells [52], a 
cell type known to be a common source of this cytokine [53], while 
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Figure 4. Macrophage transcriptional responses to Neospora infection are Ifnarl -dependent. Microarray-based expression profiling of 
Neospora (Nc) or Toxoplasma (Tg) infected wild-type or Ifnarl-/- bone marrow-derived mouse macrophages, a) heat map shows hierarchical clustering 
analysis of 833 genes differentially regulated by either Nc or Tg on either cell type relative to uninfected cells by >2-fold (FDR<5%). Two clusters of 
co-regulated genes are indicated. Each column in heatmap represents an average of duplicate arrays, and color pattern represents row Z-score. b) 
Venn diagrams showing number of differentially regulated genes for each condition, c and e) GO enrichment results for cluster 1 and 2, respectively, 
d and f) heatmaps showing top ten most strongly induced genes in cluster 1 and 2, respectively. Gene symbol and fold change of infected WT relative 
to uninfected WT control are shown for each gene. 
doi:1 0.1 371 /journal.pone.0088398.g004 



splenic red pulp macrophages produce type I interferon when 
exposed to Plasmodium-m£ected erythrocytes [54-56]. The parasite 
ligand was not identified in these studies, but the response was 
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Figure 5. The Type I interferon response to Neospora infection 
is dependent on Tlr3 and Trif. QPCR analysis of Mx1 expression in 
Neospora infected bone marrow macrophages derived from wild-type 
(WT), Myd88~'~, Myd88~ y ~ /Trir / ~ , Tlr2~ y ~/rir4^ and Tlr3~ y ~ mice. 
Error bars indicate standard deviations for three biological replicates; 
* = P<0.01. Experiment was repeated three times with similar results. 
doi:1 0.1 371 /journal.pone.0088398.g005 



found to be TW7-dependent (Toxoplasma) and Tlr9-/ Myd88- 
dependent (Plasmodium). More recendy, a type I interferon 
transcriptional signature has been noted in hepatocytes infected 
with live (but not dead) Plasmodium sporozoites, involving the 
cytosolic RNA sensor MDA5 and the MAVS adaptor protein 
[57]. Although clearly distinct from the Tlr3- and 7n/-dependent 
response to Neospora (and heat-killed Toxoplasma) reported above, 
these observations support the concept that protozoan parasites 
can be robust activators of type I interferon signaling (Fig. 4), and 
that mechanisms of innate sensing may vary depending on tissue 
or target cell type, parasite developmental stage, or other factors. 

Although we found that both parasite species are capable of 
activating type I interferon responses when killed and added to 
cells directly, the level of induction was a log order of magnitude 
higher with Neospora (Fig. 2b). This may explain why transfection of 
macrophages with Neospora RNA, but not Toxoplasma RNA, was 
sufficient to elicit a TLR3-dependent induction of type I interferon 
responsive genes. These findings suggest that there are quantitative 
and/or qualitative differences in TLR3 ligand present in RNA 
from these parasite species. Both Toxoplasma and Neospora establish 
an intracellular 'parasitophorous vacuole' (PV) distinct from the 
endo-lysosomal pathway at the time of invasion [58,59], raising 
the question of how parasite antigens, including nucleic acids, gain 
access to endosomal TLR3 or other intracellular pattern 
recognition receptors. One possibility is that TLR3 activation 
may be attributable to the release of RNA from parasites that are 
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Figure 6. Neospora RNA elicits a TLR3-dependent type I interferon response, a) QPCR analysis of Irf7 gene expression in wild-type mouse 
bone marrow-derived macrophages pretreated for 1 hr with either DMSO (control) or 100 nM bafilomycin A1 prior to 24 hr infection with live (Nc) or 
heat-killed Neospora (HK Nc). b) Irf7 expression in Myd88^ / ^ macrophages incubated with 1 ug/ml of either poly(l:C) or total RNA from HFF cells 
(host), or Neospora parasites, either alone (white) or in complex with DOTAP transfection reagent (striped), to enhance targeting of RNA to 
endosomes. c) Irf7 expression in Myd88^ / ^ macrophages incubated with 1 ug/ml of either Toxoplasma or Neospora RNA with DOTAP in cells 
pretreated with either DMSO or a small molecule inhibitor of TLR3/dsRNA complex formation. Error bars indicate standard deviations for three 
biological replicates; * = P<0.01. Experiments were repeated two or three times with similar results. 
doi:1 0.1 371 /journal. pone.0088398.g006 



phagocytosed [60] or that die after invasion, rather than those 
engaged in productive intracellular replication. Alternatively, since 
the PV is intimately associated with the host mitochondria and 
endoplasmic reticulum [6 1] , it is possible that RNA within the PV 
may access the host endo-lysosomal system (Fig. 6). PV antigens 
are known to enter the host endoplasmic reticulum for processing 
and presentation to T cells on MHC class I [62-64]. Endosomal 
vesicles may also be recruited to the PV [65,66], facilitating the 
acquisition of nutrients, and perhaps playing a role in escape from 
the host cell [67,68]. 

While type I and II interferons utilize different receptors and 
transcription factors, there is significant overlap in the genes that 
they induce [54,69], suggesting some degree of functional 
redundancy. We observed that Neospora infection drives Ifnar- 
dependent expression of genes traditionally considered to be Ifngr- 
dependent (Fig. 4), including the chemokines Cxcl9 and CxcllO, 
and the GTPases Igtp and Gbp2 (Fig. 4f), all of which have been 
shown to be important mediators of resistance to Toxoplasma 
infection [70—73]. Interestingly, Kim et ai, found that Ifnar 1 
deficiency only impacted Plasmodium infection when mice were also 
lacking the receptor for interferon-y [54], consistent with the 
notion that the strong IFN-y response typically seen in protozoan 
infections may mask a role for type I interferons. In vivo depletions 
of IFN-y were carried out to test whether IFN-y could compensate 
for loss of IFN-Ot signaling in jVfo.s/iora-infected Ifnar 1 mice (Fig. 
S3). Although WT mice treated with control IgG had no parasites 
and limited neutrophilia evident in cytospins prepared from 
peritoneal exudate cells at 7 days post-infection (Fig. S3a), WT 



mice depleted of IFN-y had abundant parasites and high 
neutrophil numbers evident (Fig. S3b), and acute mortality 
occurred at 8 dpi (data not shown). Parasites were not observed 
in Ifnar / mice treated with control IgG (Fig. S3c), however like 
their WT counterparts, all Ifnar 7 mice showed high parasitemia, 
high neutrophil infiltration and also succumbed to infection at 8 
days post-infection (Fig. S3d). Taken together, these results are 
consistent with previous reports that IFN-y is essential for control 
of Neospora [3,74], and extend these findings by showing that type I 
interferon signaling is dispensable for the control of acute high- 
dose Neospora infection. Similar results have emerged for Plasmo- 
dium infection, where parasite burden and disease pathogenesis are 
either unchanged or minimally affected by loss of type I interferon 
signaling [54,55]. 

The observation that live Toxoplasma potently suppresses the 
type I interferon pathway, and that this suppression is a dominant 
phenotype in co-infection assays with Neospora (Fig. 2) suggests at 
least two possible evolutionary scenarios: 1) that the common 
ancestor of these parasites could suppress induction of type I 
interferon and Neospora lost this ability after divergence; 2) or that 
the common ancestor could not suppress and Toxoplasma gained 
this ability after divergence. Such distinct host phenotypes could 
provide the basis for genetic screens or comparative genomic 
studies to identify the parasite factor(s) responsible for modulating 
the pathway. The fact that suppression is mediated by a soluble 
factor in the supernatant of Toxoplasma infected cells (Fig. 3) also 
raises the possibility that standard biochemical approaches could 
be used to identify and characterize a novel host or parasite factor 
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that modifies this important signaling pathway. Our observation 
that Toxoplasma suppresses TLR3-dependent type I interferon 
induction in vitro contrasts with recent findings by Koblansky et al. 
[51] in which Toxoplasma profilin elicited low levels of IFN-a from 
dendritic cells in a 77ri2-dependent manner. Given that TLR12 
signaling is MYD88-dependent, while TLR3 signaling in our 
model is not (Fig. 5), it stands to reason that Toxoplasma suppression 
may exhibit pathway specificity. In addition, Koblansky et al. 
attribute IFN-a production to plasmacytoid dendritic cells, and we 
did not test the induction or suppression of TLR3-dependent 
responses by either Toxoplasma or Neospora in this cell type. 
Interestingly, even low levels of IFN-a induced during Toxoplas- 
ma infection were critical in driving IFN-y production by NK cells 
to limit parasite replication in vivo [51], providing a possible 
evolutionary pressure for Toxoplasma parasites to limit type I 
interferon induction. Taken together with our data, a broader 
picture of innate recognition of coccidian parasites begins to 
emerge in which both protein and nucleic acid ligands are capable 
of being recognized by distinct receptors, across various cell types 
and host species. A more complete understanding of how these 
pathogens are recognized by, and interfere with, innate signaling 
pathways will help to shed light on important aspects of the host- 
pathogen relationship during infections with protozoa. 

Materials and Methods 

Ethics statement 

This study was carried out in accordance with the recommen- 
dations in the Guide for the Care and Use of Laboratory Animals 
of the National Institutes of Health. Protocols were approved by 
the Institutional Animal Care and Use (IACUC) committee of the 
University of Pennsylvania (animal welfare assurance number 
A3079-01). The University of Pennsylvania Animal Care and Use 
Programs are fully accredited by the Association for Assessment 
and Accreditation of Laboratory Animal Care International 
(AAALAC). 

Parasites, cells and mice 

RH, GT1, Prugniaud (Pru), VEG strain T. gondii; and Ncl 
[75], Nc2 [76] and NcLiv strain N. caninum were maintained by 
serial passage in human foreskin fibroblast (HFF) monolayers as 
described previously [77]. RH stably expressing mCherry were 
obtained from Dr. Anita Koshy [78]. Primary bovine fibroblasts 
were provided by Dr. Kenneth John McLaughlin (Penn). 
Cryopreserved bone marrow cells recovered from Myd88 ' , 
Myd88/Trir /_ , and Tlr2/Tlr4 _/ ~ mice were differentiated into 
macrophages as described previously [79], and infected with a 1:1 
MOI of parasites for QPCR and microarray experiments at 16— 
22 hours post-infection. Ifharl -7- and Tlr3 _/ ~ mice were 
provided by Drs. Hao Shen and Yongwon Choi (Penn), 
respectively. For mouse infections, WT and Ifnarl ' mice 
received either control IgG or anti-IFN-y (clone XMG-1.2, lmg 
LP., BioXCell) on day 0 and again at 4 days post-infection. Mice 
were infected with 3 x 10 6 NcLiv strain Neospora caninum by LP. 
injection on day 0. All mice were maintained at the University of 
Pennsylvania in accordance with Institutional Animal Care and 
Use Committee guidelines. For cytospins, 200,000 cells were spun 
onto glass slides using a Shandon Cytospin and stained with Diff- 
Quik, dehydrated, mounted with Permount (Sigma Aldrich), and 
images captured on a Nikon E600 upright microscope. 



Microarray-based expression profiling and analysis of 
Gene Ontology enrichment 

For whole genome expression microarray, total RNA was 
isolated using the RNeasy Plus kit (Qiagen). Biotin-labeled 
complementary RNA (cRNA) was generated using the Illumina 
TotalPrep RNA amplification kit. RNA and cRNA quality were 
assessed on a BioAnalyzer (Agilent). Illumina HumanRef-8 version 
3 (for HFF cells) or MouseWG-6 version 2 (for bone marrow 
macrophages) expression beadchips were hybridized with cRNA 
from two to three biological replicates per condition according to 
the manufacturer's instructions, and scanned on a beadscan unit. 
Scanned images were converted to raw expression values using 
GenomeStudio vl.8 software (Illumina). Data analysis was carried 
out using the statistical computing environment, R (v3.0.2), the 
Bioconductor suite of packages for R, and RStudio (v0.97). Raw 
data was background subtracted, variance stabilized and normal- 
ized by robust spline normalization using the Lumi package [80] . 
Differentially expressed genes were identified by linear modeling 
and Bayesian statistics using the Limma package [81,82]. Probes 
sets that were differentially regulated >2fold (FDR<5%; after 
controlling for multiple testing using the Benjamini-Hochberg 
method [83,84]) were used for hierarchical clustering and heatmap 
generation in R. Multiple probe sets per gene were collapsed to a 
single non-redundant list of expression data using the "MaxMean" 
method of the CollapseRows function in R [85]. For each gene, a 
mean expression value was calculated for biological replicates and 
used for heatmap generation (Fig. 2 and Fig. 4). Clusters of co- 
regulated genes were identified by pearson correlation using the 
hclust function of the stats package in R. Data has been deposited 
on the Gene Expression Omnibus (GEO) database for public 
access (GSE45632 and GSE45633). Gene Ontology enrichment 
analysis was conducted using either the Function Annotation 
Chart (Fig. lc) or Functional Annotation Clustering (Fig. 4) tools 
using only GO BP 'fat' terms in the Database for Visualization and 
Integrative Discovery (DAVID) [86,87]. DAVID enrichment 
scores >1.3 are equivalent to a P value<0.05 (Fig. 4c and 4e). 

Quantitative PCR (QPCR) 

cDNA was generated from total RNA using the High Capacity 
cDNA Reverse Transcription Kit (Applied Biosystems). QPCR 
was carried out using SYBR green dye and gene specific primers 
for human or mouse MX1 and IRF7. Relative transcript 
abundance was determined using the AACt method [88], using 
a standard curve of cDNA amplified with primers specific for the 
house-keeping gene GAPDH. 

Viral interference assays 

For Fig. Id, confluent HFF cells were infected with Toxoplasma 
or Neospora for 16 hr, followed by challenge with GFP-tagged 
vesicular stomatitis virus (VSV; MOI = 20 pfu/host cell) for 8 hr 
prior to assaying by fluorescence microscopy. For Fig. 3, HFF cells 
were pretreated for 4 hours with either fresh media or with 
conditioned media recovered from Toxoplasma infected HFF cells 
at 24 hr post-infection and filtered to remove parasites and host 
cell debris. Cells were then infected with Neospora for 24 hours and 
induction of IRF7 and HERC5 transcript levels was measured by 
QPCR (Fig. 3a). Supernatants from these cultures were also 
recovered, filtered again through a 0.22 |Jm filter to remove 
parasites and host cell debris, and transferred to fresh HFF cell 
monolayers prior to VSV-GFP challenge (Fig. 3b). 
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Bafilomycin treatment, DOTAP transfections and TLR3/ 
dsRNA antagonist 

Wild-type bone marrow-derived macrophages were treated with 
100 nM Bafilomycin Al for 1 hr at 37°C to block endosomal 
fusion prior to infection, RNA harvest and QPCR as described 
above. For transfection of macrophages with RNA, parasite or 
host (HFF) total RNA was isolated using the miRNeasy kit 
(Qiagen). Myd8S~'~ bone marrow-derived macrophages were 
used in order to minimize induction of innate signaling pathways 
other than Tlr3. 1 u,g total RNA or poly(LC) (Imgenex) was added 
to each chamber of a 24-well plate containing 1 ml of media and 
1x10 s Mjd88 1 macrophages. For more efficient targeting to 
endosomes, 1 (Xg of parasite RNA, host RNA, or poly(LG) was 
mixed with 10 u.g l,2-dioleoyl-3-trimethylammonium-propane 
(DOTAP) liposomal transfection reagent (Roche Applied Science). 
Cells were incubated for 18 hr with transfection mixture and RNA 
was isolated for QPCR analysis as described above. To test a role 
for TLR3 in recognition of parasite RNA, a thiophenecarbox- 
amidopropionate small molecule inhibitor of the TLR3/dsRNA 
complex (Calibochem) was added to cultures at 50 nM for lhr at 
37°C before transfection of cells with RNA/DOTAP. 

Assaying for presence of parasite RNA virus 

To test for presence of an endogenous RNA virus in the 
parasite, total RNA was extracted from purified Toxoplasma and 
Neospora tachyzoites using Trizol. Nuclease digestions were carried 
out using 60 u.g of total RNA in 1 5 u.1 of volume containing 3 
units/Ml of SI nuclease (Invitrogen), alone or following Dnase 
treatment and RNA purification, and incubated at 37°C for 45 
minutes. Samples were analyzed on 0.8% native agarose in 1 x 
TAE buffer at 4°C. SI nuclease-treated Leishmania guyanensis (strain 
M4147), which contains the endogenous RNA virus LRV1, was 
used as a positive control. 

Fluorescence microscopy to assess parasite 
internalization 

HFF cells grown to confluency on coverslips were infected with 
a 3:1 MOI of either live or heat-killed (56°C for 30 minutes) RH 
strain Toxoplasma expressing mCherry [78]. 18 hours after 
infection, coverslips incubated with 50 nM LysoTracker (Molec- 
ular Probes) for 10 minutes at 37°C were live mounted, and 
immediately imaged using a Leica DMI400 with Yokagawa 
spinning disk confocal system. Images were overlayed and 
analyzed using Image J software, vl.47. 

Statistical Analysis and data visualization 

All experiments were repeated 2-4 times, and means and 
standard deviations were calculated from biological replicates. 
Significance was determined using a Student's t-test. Statistical 
analysis was carried out using Prism 4 (GraphPad Software) and 
data visualization using DataGraph 3.0 (Visual Data Tools). 

Supporting Information 

Figure SI Assaying for double stranded RNA in Toxo- 
plasma and Neospora. 0.8% agarose gel electrophoresis and 
ethidium bromide staining of Leishmania guyanensis (Lg), Toxoplasma 
gondii (Tg) and Neospora caninum (Nc) total RNAs treated with A) SI 
nuclease, B) DNase I and SI nuclease, C) untreated total RNA of 
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QPCR analysis of the expression of the antiviral gene Mxl in (a) 
HFF cells infected with Neospora tachyzoites at various multiplicity 
of infection (MOI), and (b) bovine primary fibroblasts infected with 
either Toxoplasma (Tg) or Neospora (Nc). Similar results were 
obtained for Irf7 (not shown). Error bars indicate standard 
deviations for two biological replicates; * = PS0.01. 
(TIFF) 

Figure S3 IFN-y, but not type I interferons, are essential 
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Toxoplasma strains. 66 jV"««/wra-specific antiviral genes (Figure lb) 
are shown on separate tab of this Excel file. 
(XLSX) 

Table S2 Illumina microarray average log2 expression 
values from Figure 4a showing 833 mouse genes 
differentially regulated >2-fold (FDR<5%) following 
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